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Abstract 

We consider phase coexistence and criticality in a thin-film Ising magnet 
with opposing surface fields and non-planar (corrugated) walls. We show 
that the loss of translational invariance has a strong and unexpected non- 
linear influence on the interface structure and phase diagram. We identify 4 
non-thermodynamic singularities where there is a qualitative change in the 
interface shape. In addition, we establish that at the finite-size critical point, 
the singularity in the interface shape is characterized by two distint critical 
exponents in contrast to the planar case (which is characterised by one). 
Similar effects should be observed for prewetting at a corrugated substrate. 
Analogy is made with the behaviour of a non-linear forced oscillator showing 
chaotic dynamics. 

PACS numbers: 68.45. Gd, 68.10.-m, 68.35. Rh 
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There are a number of well studied examples of fluid interfacial phenomena for planar 
systems in which surface phases with distinct adsorptions co-exist along a line of first-order 
phase transitions which terminates at a surface critical point. Examples include the pre- 
wetting transition associated with first-order wetting [|l[] and also interfacial localization 
in thin-film magnets (with opposite surface fields) associated with confinement effects at 
critical wetting 0. In both cases, the difference in adsorption between the two phases 
vanishes continuously as the critical point, signifying the end of two-phase coexistence, is 
approached. This second order phase transition is characterized by the critical exponents 
belonging to the two dimensional Ising universality class (for three dimensional bulk systems) 
since the adsorption difference acts as a scalar order parameter In this letter, we 

describe a wealth of new interfacial structural changes and singularities which emerge when 
the analogous phenomena are considered in slightly non-planar geometries and which are 
intimately associated with non-linear behaviour. In addition to a shift in the finite-size 
(FS) critical point (compared to the planar confined system), the shape of the non-planar 
interface undergoes a number of structural changes as we move along and beyond the line 
of coexistence. This behaviour has no counterpart in the planar geometry, and has not 
been reported previously. Moreover, as the shifted surface critical point is approached, the 
function describing the shape of the non-planar interface shows non-analyticities which are 
characterised by two critical exponents. Whilst one of these appears to be identical with the 
usual critical exponent describing the singularity in the total (or average) adsorption, the 
general identification of the second is a more difficult problem, although scaling arguments 
(consistent with our explicit results) suggest that its value is related to the energy density. 
Our predictions are based on a detailed numerical analysis of a simple mean-field (MF) model 
of interfacial behaviour which we believe is qualitatively correct beyond MF approximation 
(in three dimensions). These are rather dramatic effects emanating from the introduction 
of a slight non-planar perturbation to the interface can be viewed profitably by making 
analogy with the classical mechanics of an extremely sensitive non-linear dynamical system 
exhibiting chaotic behaviour ||. As we will see, the interface behaviour may be elegantly 
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portrayed as the temperature evolution of a phase plane plot, similar to that employed in 
dynamical systems, allowing us to distinguish different qualitative types of interface shape 
separated by non-thermodynamic singularities. 

To begin, we recall the relevant properties and phase diagram of the planar system prior 
to a discussion of the non-planar generalization. The transition that we concentrate on 
occurs in a thin-film magnet with opposing surface fields but the phenomena are generic 
to other situations such as prewetting at a planar wall. Consider then an Ising-like thin 
film magnet of width L z and infinite transverse area in zero bulk field and below the bulk 
critical temperature Tf VLK with surface fields h\ and hi = —hi acing on the spins in the 
z = and z = L z planes respectively. We further suppose (through a judicial choice of 
surface enhancement H) that in the semi-infinite limit L z — > oo each surface undergoes a 
critical wetting transition at temperature T w . For such a system, MF and simulation 
studies || show that the finite size phase diagram is dominated by wetting effects which 
are able to suppress bulk-like coexistence over a large temperature regime. At sufficiently 
low temperatures T <T C (L Z ), with the finite-size critical temperature satisfying T C (L Z ) <T W , 
phase coexistence is possible between phases corresponding to an interface being bound 
to either wall. As the temperature increases, the interface position moves continuously to 
the middle of the system and for T > T C (L Z ) only one phase is possible. Thus, in the 
temperature window T^ ULK >T>T W , the FS effects suppress bulk-like phase coexistence 
for all L z . This temperature range is also characterised by a near soft-mode phase since 
the transverse correlation length £11 is extremely large due to capillary- wave like excitations. 
These features can be most easily understood using a simple effective interfacial Hamiltonian 
model @: 



where £(r) is the collective co-ordinate describing the interface position at vector displace- 
ment r = (x, y) along the wall and £ is the stiffness coefficient of the up-spin-down-spin 
interface. The total finite-size binding potential W(£; L z ) acting on the interface (whose 




(1) 
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minima determine the MF location/s of the interface) is the sum of the two contributions 
from each wall: 



W(£; L z ) = WM + W^{L Z - £) (2) 

where W / oc (£) is the appropriate semi-infinite binding potential for the ranges of forces in 
the model. For systems with short ranged forces this is usually specified as [ij 

W(£) = a (T - T w ) e~ Ki + b e~ 2Ki ; £ > (3) 

with a , b a positive constants and k being the inverse bulk correlation length. For T < T C (L Z ) , 
with T C (L Z ) = T W — A{b / a )e~ KL ^ 2 in MF approximation, the total potential W(£; L z ) has a 
double well structure with two equal minima at £ n < L z /2 and = L z —£ n . As T -^T C (L Z )~ , 
the adsorption difference Ar = Am (L z /2 — £ n ) (with m the bulk magnetization) vanishes 
like Ar ~ (T C (L Z ) — TJ 1 / 2 , corresponding to a standard order- disorder transition. For 
T>T C (L Z ), the potential W(£;L Z ) has only one minimum at £ n = L z /2 and the correlation 
length £|| ~ e K W 4 . Interestingly, most of these quantitative MF predictions are confirmed by 
extensive Monte Carlo simulation studies which established that the true asymptotic critical 
regime where we can expect Onsager-like behaviour Ar ~ (T C (L Z ) — T c ) l / S is extremely small 
0]. All these facts support MF theory as an excellent quantitative description of the thin 
film system. 

We now wish to consider the MF phase diagram for the analogous phase transition in a 
slightly non-planar geometry. We will take as our starting point the simplest possible phe- 
nomenological model of this system which generalises ([!]) and suppose that the configuration 
energy is specified by 



H[£;z^\z^]= Jdr 



^(V£f + W(£;L z ,zV,z^) 



(4) 



where z!£'(r) and zffl{r) describe the (small) deviations of the walls near the z = and 
z = L z respectively and W(£; L z , z$,z$) =W OQ (£- z$) + W^L, + z^ - £). Whilst the 
model could certainly be improved by including further coupling terms involving V£ ■ Vz w 



with associated position dependent (stiffness) coefficients, we do not expect these to make 
any significant difference to the interfacial behaviour described here ||. In any case, even 
with the further assumption of corrugated walls such that z^ and z$ only depend on 
a single-coordinate (x say), the interfacial behaviour generated is sufficiently complex to 
warrant attention within the simple model above. Writing z!£'(x) = ay2sm(qx), we have 
considered the geometry for which zffl = zffl although, of course, many other choices are 
possible 0. The r.m.s. width a and wavelength L x = 27r/q of the wall corrugation are 
assumed to be small and large respectively in comparison with the bulk correlation length. 
With these assumptions, the wetting transition remains second order and located at T w in 
the semi-infinite limit [0. 

The equilibrium non-planar interfacial profile/s t v (x) satisfies the Euler-Lagrange equa- 
tion 

£ i v {x) = W^(£ v - z«) - W^(L Z + - 4) (5) 

where dot and prime signify differentiation w.r.t. x and argument respectively. Periodic 
boundary conditions are imposed after a large multiple of wavelengths L x . Two preliminary 
remarks are as follows: firstly, the Euler-Lagrange equation is inversion symmetric so that if 
£ u (x) is a solution, £*(x) =L Z —£ v (x+ir/ q) is also a solution with the same free energy and is 
distinct from t v (x) in the two-phase regime. Secondly, we have established numerically that 
the stable phases all have the same wavelength as the wall corrugation L x . However, this 
is not the case for the metastable states ||. Finally, we note that an elegant description of 
the interfacial shape is afforded by a reduce phase plane plot £/y2a vs. £/\/2aq and helps 
distinguish different types of structural regimes. A section of the equilibrium phase diagram, 
with suitable reduced units |§, is shown in Fig. 1 and shows a critical line (corresponding 
to an order-disorder transition) and four non-thermodynamic singularities where there is 
a qualitative change of interfacial structure. In this way, we are able to distinguish five 
different interfacial types (see Fig. 2). 

Phase coexistence and order is most easily revealed through the mean interfacial height 
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£ = J- / L rfx (6) 
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L x Jo 

which is single valued (£ Q = L z /2) in the disordered regime above the critical temperature 
T C (L Z) a, q), but is double valued (with £* = L z — £ ) in the order regime, analogous to £ n and 
£* for the planar system. Our numerics indicate that the singularity in £ Q is of the expected 
type: 

if t > 

(7) 

if t<0 

where we have introduced the scaled temperature variable t = (T c (L z , a, q) — T)/T. In addi- 
tion to the mean interface height, however, the shape function shows a number of qualitative 
changes with temperature. At very low temperatures, the interface is closely bound to one 
of the walls and follows the corrugation (See Fig. 2(A)). Over one period L x , the graph 
£ v {x) has one maximum and one minimum which are in phase with the wall function zffl(x). 
For this case, the phase plane plot is a simple loop. Nevertheless, notice that its form is not 
precisely circular, indicating that non-linear effects are important even when the interface is 
close to the wall. On increasing the temperature, the interface smoothly deforms and shows 
a number of non-thermodynamic singularities where the minima and maxima of the graph 
undergo a series of bifurcations. These reveal themselves as the appearance/disappearance 
of loops in the phase plane portrait as illustrated in Fig. 2(B) which also shows the locus 
of the maxima/minima with temperature (Fig. 2(C)). Corresponding profiles are shown in 
Fig. 2(A). Two counter intuitive features are worth emphasising here. Firstly, there are two 
regimes, II and IV, where the interface shape has two and three maxima per wavelength 
of the wall corrugation. Secondly, in the vicinity of the order-disorder transition, regime 
III, the interface shape is similar to the wall (i.e. there is only one max/min pair per pe- 
riod) but is out of phase with it. Finally, at high temperatures above the two super critical 
non-thermodynamic singularities, the interface shape returns to that of a simple sinusoidal- 
like function in phase with the wall and the phase portrait is basically a circle of radius 
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(1 + q 2 ^y) _1 centred at L z /2. 

Next, we focus on the singularity in the shape profile at the order-disorder transition. 
We have established that the stable phase/s can be represented by a Fourier series 

£ v (x)=£ + a 1 wn(qx)+ (8) 

oo 

{ a 2k+i sin ((2k + l)qx) + 7 2fc cos(2/cqx)} 

fc=i 

throughout the phase diagram. In this expression i is the mean interface position (given by 
Eq. (§)), whilst the second term is the harmonic response to the wall corrugation. The final 
term represents the higher order harmonic excitations arising from the non-linearity of the 
Euler-Lagrange equation and are responsible for the complicated evolution of the interface 
structure with temperature. We stress that, without this term (i.e. just considering linear 
response), the phase plane portrait would be simply circular. Note that there are not 
even sine terms and no odd cosine terms. The temperature dependence of the two sets of 
coefficients {<J2k-i} and {72k} is extremely involved but near T c (L z ,a, q) only two types of 
singularity are observed in our numerical analysis (See Fig. 3). The coefficients {72^} all 
vanish above T C (L Z , a, q) and behave precisely as the mean order-parameter 4* — £ Q , i.e. they 
are characterised by the usual MF order-parameter critical exponent /? = 1/2. In contrast, 
the terms {a2k-i} all have a cusp-like singularity 

ffa-i-ffLi- \t\°; i^0 ± (9) 

where <J2k-i is the value at criticality and the critical exponent 6 = 1. There is no analogy 
of this singularity in the planar system. Furthermore, whilst it is natural to identify the 
cosine term singularities with the order-parameter exponent (3 of the d — 1 dimensional bulk 
universality class ((3 = 1/2 in MF, (3 = 1/8 beyond MF for three dimensional thin films), 
a similar identification for 6 is not so obvious. Nevertheless, we have constructed scaling 
arguments which suggest that 6 = 1 — a, where a is the specific heat critical exponent, 
consistent with our numerical results ||. Similarly, we have also established that, for fixed 
L z , the critical line is consistent with the scaling law 
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T c (L z ,a,q) - T C (L Z , 0, 0) 




(10) 



where A is an appropriate scaling function. This behaviour can be understood using finite- 
size scaling ideas with MF critical exponents and indicates that the effective width of the 
system is reduced by corrugation ||. 

To finish our article, we make some pertinent remarks. Firstly, the interfacial structural 
changes reported here are not peculiar to short-ranged forces with the exponential binding 
potential Eq. @, and also emerge if long-ranged forces are considered instead |§. Also, we 
emphasise that in previous studies of the effect of roughness on wetting transitions most 
authors have considered binding potentials with a single minimum which do not exhibit 
the same subtle non-linear behaviour discussed here J?H§. Next, we note that on making 
a change of variable r]{x) =£(x) — L z —zffl and expanding to appropriate (cubic) order, the 
Euler-Lagrange equation can be written 



where t oc (T C (L Z ) — T) and u is positive in the region of interest. This is essentially 
equivalent to the Duffing equation of a soft-polynomial oscillator (without a damping term) 
which is known to yield extremely rich (including chaotic) dynamics ||. In this context, the 
non-thermodynamic singularities described above are analogous to the harmonic excitations 
of the non-linear oscillator (however, this analogy does not shed any light on the nature of 
the singularities near the order-disorder transition and their identification beyond MF). 

In summary, we have shown from a simple MF model of interfacial behaviour in a slightly 
non-planar geometry that new types of structural phase changes and additional critical 
singularities can emerge which are intimately related to non-linear phenomena. Similar 
behaviour is also expected for pre- wetting at a non-planar substrate ||. Also of interest is 
the structure of metastable states in this system which we do not discuss here |§. We believe 
that future studies of improved models which include thermal fluctuations and different type 
of non-planarity will also uncover new structural and fluctuation related behaviour. 



E f] = —trj + urf + aq 2 sin(qx) 



(11) 
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FIGURES 



FIG. 1. Phase diagram for L z = 10 and q = 2n/W in reduced units ||. The solid line 
separates the ordered and the disordered phases. The dashed lines show the location of the 
non-thermodynamic singularities and divide the phase diagram into 5 regions. 



FIG. 2. Behaviour of the system for L z = 10, a = 1.5 and q = 2n/10 showing the shape of the 
interface in the different regimes (A) and their corresponding phase portraits £ vs. £ (B). The circle 
represents the point x = 0. For clarity, scales related to £ are ommited but can be checked in Fig. 3. 
The loci of the interface minima and maxima are represented as a funcion of the temperature (C). 
The FS critical temperature T c (L z ,a, q) ~ 0.845 is represented by a thin line and is within regime 

in. 



FIG. 3. Behaviour of the coefficients L z /2 — £ Q , ci, 72 and <r 3 of Eq. || near T C (L Z , o, q) for 
L z = 10, o=1.5 q i = 27r/10. They are multiplied by 1, WL X , 10 2 L X and lO 3 !^, respectively. 
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